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Abstract We give an overview of the magnetic catalysis phenomenon. In the frame- 
work of quantum field theory, magnetic catalysis is broadly defined as an enhance- 
ment of dynamical symmetry breaking by an external magnetic field. We start from 
a brief discussion of spontaneous symmetry breaking and the role of a magnetic 
field in its a dynamics. This is followed by a detailed presentation of the essential 
features of the phenomenon. In particular, we emphasize that the dimensional re- 
duction plays a profound role in the pairing dynamics in a magnetic field. Using the 
general nature of underlying physics and its robustness with respect to interaction 
types and model content, we argue that magnetic catalysis is a universal and model- 
independent phenomenon. In support of this claim, we show how magnetic catalysis 
is realized in various models with short-range and long-range interactions. We argue 
that the general nature of the phenomenon implies a wide range of potential applica- 
tions: from certain types of solid state systems to models in cosmology, particle and 
nuclear physics. We finish the review with general remarks about magnetic catalysis 
and an outlook for future research. 



1 Introduction 



The magnetic catalysis is broadly defined as an enhancement of dynamical symme- 
try breaking by an external magnetic field. In this review, we discuss the underlying 
physics behind magnetic catalysis and some of its most prominent applications. 
Considering that the ideas of symmetry breaking take the center stage position in 
many branches of modern physics, we hope that this review will be of interest to a 
rather wide audience. 
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In particle and nuclear physics, spontaneous symmetry breaking is commonly 
used in order to explain the dynamical origin of the mass of elementary particles. 
In this context, the idea was realized for the first time over 50 years ago by Nambu 
and Jona-Lasinio 11 1491 11501 . who suggested that "the nucleon mass arises largely 
as a self-energy of some primary fermion field through the same mechanism as the 
appearance of energy gap in the theory of superconductivity." As we now know, 
the analogy with superconductivity is not very close and the description of chiral 
symmetry breaking in terms of quarks may be more natural than in terms of nucle- 
ons. However, the essence of the dynamical mass generation was captured correctly 
in Refs. 1149111501 . In fact, with the current state of knowledge, we attribute most 
of the mass of visible matter in the Universe to precisely this mechanism of mass 
generation, which is associated with breaking of the (approximate) chiral symmetry. 

The conceptual knowledge that the mass can have a dynamical origin opens myr- 
iads of theoretical possibilities that would appear meaningless in classical physics. 
For example, keeping in view the above mentioned mechanism of mass generation 
through chiral symmetry breaking, it is reasonable to suggest that the masses of 
certain particles can be modified or even tuned by proper adjustments of physical 
parameters and/or external conditions. 

One of the obvious knobs to control the value of the dynamical mass is an exter- 
nal magnetic field. In addition to be a good theoretical tool, magnetic fields are also 
relevant to many applications. For example, they are commonly present and play an 
important role in such physical systems as the Early Universe 11811 144, 30, 181 I951 . 
heavy ion collisions 1122111771 . neutron stars (1801 [381 , and quasi-relativistic con- 
densed matter systems like graphene (11531118 8l . 

As we discuss in detail in this review, the magnetic field has a strong tendency to 
enhance (or "catalyze") spin-zero fermion-antifermion condensates. Such conden- 
sates are commonly associated with breaking of global symmetries (e.g., such as 
the chiral symmetry in particle physics and the spin-valley symmetry in graphene) 
and lead to a dynamical generation of masses (energy gaps) in the (quasi-)particle 
spectra. The corresponding mechanism is called magnetic catalysis l99l . 

It should be emphasized that, in a striking contrast to its role in superconduc- 
tivity, the magnetic field helps to strengthen the chiral condensate. There are many 
underlying reasons for its different role. Unlike the superconductors, the ground 
state with a nonzero chiral condensate shows no Meissner effect. This is because the 
chiral condensate can be thought of as a condensate of neutral fermion-antifermion 
pairs, not charged Cooper pairs that can give rise to supercurrents and perfect dia- 
magnetism. Also, in a usual Cooper pair, the two electrons have opposite spins and, 
therefore, opposite magnetic moments. When placed in a magnetic field, only one 
of the magnetic moments can minimize its energy by orienting along the direction 
of the field. The other magnetic moment will be stuck in a frustrated position point- 
ing in the opposite direction. This produces an energy stress and tends to break the 
Cooper pair. (Note, however, that the orbital motion plays a much more important 
role in breaking nonrelativistic Cooper pairs.) In a neutral spin-zero pair, in contrast, 
the magnetic moments of the fermion (with a fixed charge and spin) and the an- 
tifermion (with the opposite charge and spin) point in the same direction. Therefore, 
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both magnetic moments can comfortably align along the direction of the magnetic 
field without producing any frustration in the pair. (Also, in relativistic systems the 
fermion-antifermion condensate is not destroyed by the orbital motion.) 

The above explanation of the role that the magnetic field plays in strengthening 
the chiral condensate is semi-rigorous at best and does not capture all the subtleties 
of the dynamics behind magnetic catalysis (e.g., completely leaving out the details 
of the orbital motion). It does demonstrates, however, how the magnetic field can 
have, at least in principle, so drastically different effects on the dynamical gener- 
ation of mass on the one hand and on superconductivity on the other. (It may be 
curious to mention here that, in cold dense quark matter, it is possible to obtain 
color superconducting states, in which diquark Cooper pairs are neutral with re- 
spect to the in-medium (but not vacuum) electro-magnetism 1831 171. In such quark 
matter, the in-medium magnetic field is not subject to the Meissner effect and, in 
fact, can enhance color superconductivity |[T52l|75ll62ll^[l39ll54ll55l[T84l.) 

The early investigations of the effects of strong magnetic fields on chiral sym- 
metry breaking in (2 + 1)- and (3 + l)-dimensional models with local four-fermion 
interactions have appeared in late 1980s and early 1990s lfT27l [T79l [1291 [1281 [1631 
1133111261 . In these studies, it was already found that a constant magnetic field sta- 
bilizes the chirally broken vacuum state. 

The explanation of the underlying physics was given in Ref. [1991 . where the 
essential role of the dimensional reduction, D — > D — 2, in the low-energy dy- 
namics of pairing fermions in a magnetic field was revealed. As a corollary, it 
was also established that the presence of a magnetic field leads to the generation 
of a dynamical mass even at the weakest attractive interaction between fermions 
ll99l 1 1 001 1 1 1 1 1 1 021 . The general nature of the underlying physics was so compelling 
that it was suggested that the corresponding dynamical generation of the chiral con- 
densate and the associated spontaneous symmetry breaking in a magnetic field are 
universal and model-independent phenomena. To emphasize this fact, the new term 
"magnetic catalysis" was coined 11991 . 

The model-independent nature of magnetic catalysis was tested in numerous (2 + 
1)- and (3 + 1) -dimensional models with local four-fermion interactions lfl3l 11301 

ss HS21 tun insi eh m usa eu eei nm nm na w\ m eh, including 

models with additional gauge interactions 01 191 . higher dimensional models 0841 . 
,jV = 1 supersymmetric models 042l . quark-meson models |8]|9), models in curved 
space lUEDHH and QED-like gauge theories H01lll36lH56lH57l[TT0l[T04l l5Tl 
[1351 H05l ED EU DSD E MB [163] H2 El ■ The realization of magnetic catalysis 
was investigated in chiral pertubation theory 11741 13211331 and in QCD 0144111211 . 
as well as in a models with the Yukawa interaction [61 60 59 . 43 ) . There are studies 
of magnetic catalysis using the methods of the renormalization group 0169lll64ll74l . 
lattice calculations |[52j |24] [28j [34j [35] [J6] [T4l [T51 and holographic dual modelss of 
large-N gauge theories £[85] M M USD SS S3 SU S3 M EEO S3 ED ED M 
SI- Similar ideas were extended to solid state systems describing high-temperature 
supersonductivity 11371 11681 11931 l57l l58l 11911 , highly oriented pyrolitic graphite 
11241 [1231 [881 . as well as monolayer (98] [87] [89] EU ETJ [172] ED and bilayer 
ll93l l92l ED |90) graphene in the regime of the quantum Hall effect. Finally, the 
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generalization of magnetic catalysis was also made to non-Abelian chromomagnetic 
fields |[T83l[l3T][T73l|97l|4Tl[39l[l90l, where the dynamics is dimensionally reduced 
by one unit of space, D — >• D — 1 . For earlier reviews on magnetic catalysis, see 
Refs. 04311961. 



2 The Essence of Magnetic Catalysis 

As already mentioned in the Introduction, the essence of magnetic catalysis is in- 
timately connected with the dimensional reduction, D — > D — 2, of charged Dirac 
fermions in the presence of a constant magnetic field. In this section, we discuss in 
detail how such a dimensional reduction appears and what implications it has for 
the spontaneous symmetry breaking. 



2.1 Dimensional reduction in a magnetic field 

Before considering a fully interacting theory and all details of the dynamics re- 
sponsible for the generation of the chiral condensate and the symmetry breaking, 
associated with it, let us start from a free Dirac theory in a constant external mag- 
netic field. This appears to be a perfect setup to understand the kinematic origin of 
the dimensional reduction, D — > D — 2. 

2.1.1 Dirac fermions in a magnetic field in 3 + 1 dimensions 

Let us start by reviewing the spectral problem for charged (3 + 1 ) -dimensional Dirac 
fermions in a constant magnetic field. We assume that the field is pointing in the 
positive x 3 -direction. The corresponding Lagrangian density reads 

^ = W{iy^ l D ll -m) x ¥, (1) 

where the covariant derivative = <9 M — j'eAff' depends on the external gauge field. 
Without loss of generality, the external field A™' is taken in the Landau gauge, A™' = 
(0, -A ext ), where 

A ext = (0,Bx\0), (2) 

and B is the magnetic field strength. By solving the Dirac equation of motion, one 
finds the following energy spectrum of fermions (3) : 

E^ +l) (p 3 ) = ±^m^ + 2\eB\n + (p 3 )\ (3) 
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where n = 0, 1,2, ... is the Landau level index. It should be noted that the Lan- 
dau level index n includes orbital and spin contributions: n = k + s + j, where 
k = 0,1,2,... is an integer quantum number associated with the orbital motion, 
while s = ±5 corresponds to the spin projection on the direction of the field. [For 
the orbital part of the wave functions, see Eq. (IHZt in the Appendix.] Considering 
that the energy depends only on n, we see that the energy of a quasiparticle in orbital 
state k and spin s = +^ is degenerate with the energy of a quasiparticle in orbital 
state k+ 1 and spin s = — A. The lowest Landau level with n = is special: it corre- 
sponds to the lowest orbital state k = and has only one spin projection s — — i. The 
letter, in particular, implies that the lowest Landau level is a spin polarized state. 

On top of the spin degeneracy of higher Landau levels (« > 0), there is an ad- 
ditional (infinite) degeneracy of each level with a fixed « and a fixed value of the 
longitudinal momentum p?,. It is connected with the momentum p2 G R, which is a 
good quantum number in the Landau gauge utilized here. As follows from the form 
of the orbital wave functions in Eq. d82i i. the value of — p2/\eB\ also determines the 
location of the center of a fermion orbit in the x'-direction. A simple analysis J5] 
shows that the area density of such states in the perpendicular x l x 2 -plane is @ for 

n = and ^p- for n > (here the double spin degeneracy of the higher Landau levels 
is accounted for). 

When the Dirac mass is much smaller than the corresponding magnetic en- 
ergy scale (i.e., m -C \f\eW\), we find that the low-energy sector of the Dirac 
theory is determined exclusively by the lowest Landau level (n = 0). As we see 
from Eq. (0, the corresponding spectrum of the low-energy excitations is given by 
E(pi) = ±\/m 2 + (p?,) 2 , which is identical to the spectrum of a (1 + 1) -dimensional 
quantum field theory with a single spatial coordinate, identified with the longitudinal 
direction. This spectrum of the low-energy theory confirms the obvious kinematic 
aspect of the dimensional reduction, 3 + 1 ^ 1 + 1, in a constant magnetic field. 

From the physics viewpoint, the dimensional reduction is the result of a partially 
restricted motion of Dirac particles in the x l x 2 -plane perpendicular to the magnetic 
field. The effect can be seen already at the classical level in the so-called cyclotron 
motion, when the Lorentz force causes charged particles to move in circular orbits 
in x l x 2 -plane, but does not constrain their motion along the x 3 -direction. A very im- 
portant new feature at the quantum level is the quantization of perpendicular orbits. 
Without such a quantization, the clean separation of the low-energy sector, domi- 
nated exclusively by the lowest Landau level, would not be possible. 

It should be noted that the spin also plays an important role in the dimensional 
reduction of Dirac particles. If the spin contribution were absent (s = 0), the energy 
of the lowest Landau level would scale like v/Jefif, which is not vanishingly small 
compared to the energy of the next Landau level ^/3|eB|. Then, a clean separation 
of the lowest Landau level into a dimensionally reduced, low-energy sector of the 
theory would become unjustified and meaningless. 
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2.1.2 Dirac fermions in a magnetic field in 2 + 1 dimensions 

It is straightforward to obtain the spectrum of charged Dirac fermions also in 
2 + 1 dimensions. The vector potential in the Landau gauge takes the form: A ext = 
(0,-Bx 1 ) . In the absence of the longitudinal direction x 3 , the magnetic field B is not 
an axial vector, but a pseudo-scalar. Concerning the Dirac algebra in 2 + 1 dimen- 
sions, there exist two inequivalent irreducible representations, given by 

7° = 03, y l = iOi, f = io 2 , (4) 

and 

7° = -ff3, Y l =-i<7u J 2 = —ioz, (5) 

where (7, are the Pauli matrices. In each of these representations, the nature of the 
lowest Landau level is somewhat unusual: it has either only a particle state (with 
a positive energy £0 = m) or only an antiparticle state (with a negative energy 
Eq = —m). Such an asymmetry in the spectrum is known to induce a Chern-Simons 
term in the gauge sector of the theory 1151111621 . In order to avoid the unnecessary 
complication, it is convenient to use the following reducible representation instead: 




(Incidentally, in the low-energy theory of graphene, which is a real quasi-relativistic 
system in 2 + 1 dimensions, such a reducible representation appears automatically 
M170I .) The corresponding Dirac spectrum reads 

+1) = ±^m 2 + 2\eB\n. (7) 

As we see, this is very similar to the (3 + 1) -dimensional result in Eq. (O, except 
for the missing dependence on the longitudinal momentum pi. 

Repeating the same arguments as in the (3 + l)-dimensional case, we find that 
the low-energy sector of the Dirac theory in 2 + 1 dimensions is also determined by 
the lowest Landau level. Here we assume again that the Dirac mass is much smaller 
than the corresponding Landau energy scale (m <C v/jeB[) in order to insure a clear 
separation of the low- and high-energy scales. 

Just like in the higher dimensional case, all Landau levels are (infinitely) degen- 
erate. In particular, the number of degenerate states per unit area is @ in the lowest 
Landau level. A special feature of the (2 + 1 ) -dimensional theory is a discrete, rather 
than continuous spectrum of excitations. In the absence of the x 3 -direction and the 
associated quantum number /?3, all positive energy states in the lowest Landau level 
have the same energy Eq = m. Moreover, when m 0, this energy goes to zero and 
becomes degenerate with the negative energy states £0 = — m. In this limit, there is 
an infinite vacuum degeneracy even if the condition of charge neutrality may favor 
a state with exactly half-filling of the lowest Landau level. It should be expected, 
however, that taking into account any type of fermion interaction will lead to a well 
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defined ground state, in which the interaction energy is minimized. One can even 
make an educated guess that the corresponding ground state should be a Mott-type 
insulator with a dynamically generated mass/gap. 

Another special and rather unusual feature of the (2+ 1) -dimensional Dirac 
fermions in a magnetic field is a spontaneous symmetry breaking, which is man- 
ifested by a nonzero "crural" condensate (WW) already in the free theory. To see 
this, let us make use of the proper-time representation of the fermion propagator in 
the magnetic field, see Eq. ( 110U in the Appendix. In the limit of a small bare mass 
(mo — > 0), we easily derive the following (regularized) expansion for the condensate: 



It can be shown that the first term, which remains nonzero even in the massless limit, 
comes from the lowest Landau level. At first sight, this may appear to be a very sur- 
prising result. Upon a closer examination, one finds that this condensate is directly 
connected with a nonzero density of states and a nonzero spin polarization in the 
lowest Landau level of the free Dirac theory. The result is unambiguous only after 
specifying the sign of the bare mass parameter, which is also typical for spontaneous 
symmetry breaking. 

In connection with the result in Eq. ([8j, it may be useful to recall that the chirality 
is not well defined in the (2+ 1) -dimensional space. However, as we will discuss 
in Sec. 12.21 the condensate (WW) is still of interest because it breaks another global 
symmetry that has a status similar to that of the conventional chiral symmetry. 



2.2 Magnetic catalysis in 2+1 dimensions 

Now, let us consider a Nambu-Jona-Lasinio (NJL) model in 2 + 1 dimensions, in 
which the magnetic catalysis of symmetry breaking is realized in its simplest possi- 
ble form Il99lll02| . When using the reducible representation of Dirac algebra, given 
by Eq. (O, one finds that the kinetic part of the massless Dirac theory is invariant 
under a global U (2) flavor symmetry. The generators of the symmetry transforma- 
tions are given by 7b =I,T\ = y 5 ,72 = jy 3 , and T3 = y'y 5 , where y 5 = iy / 1 / 2 / 3 . 
A dynamical Dirac mass will break this U (2) symmetry down to the U(l) x f/(l) 
subgroup with generators 7b and 73 . 

The NJL-type Lagrangian density, with the interaction term invariant under the 
U (2) flavor symmetry, can be written down as follows: 




(8) 



Jz? = Wiy f "D tl W+^ (WW) 2 + (WifW) 2 + (WfW) 2 , 



(9) 
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where G is a dimensionfull coupling constant. This theory is nonrenormalizable, but 
can be viewed as a low-energy effective theory with a range of validity extending 
up to a certain ultraviolet energy scale set by a physically motivated choice of the 
cutoff parameter A . 



2.2.1 Weak coupling approximation 

As usual in problems with spontaneous symmetry breaking, we use the method of 
Schwinger-Dyson (gap) equation in order to solve for the dynamical mass parame- 
ter. We assume that the structure of the (inverse) full fermion propagator is the same 
as in the free theory, but has a nonzero dynamical Dirac mass m, 

S- l (x,x') = -i [ifd t - -m] 8 3 (x-x'). (10) 

In the mean-field approximation, the dynamical mass parameter satisfies the follow- 
ing gap equation: 

m = Gtr[S(x,x)]. (11) 

To leading order in weak coupling, G — > 0, this equation can be solved perturba- 
tively. Indeed, by substituting the condensate calculated at niQ — > in the free theory, 
see Eq. (O, into the right-hand side of the gap equation ( fTTT i. we obtain 

l«NGJgi. (12) 

As we see, the dynamical mass is induced at any nonzero attractive coupling. 

The massless NJL model in Eq. (0 is invariant under the flavor U (2) symmetry. 
The generation of a Dirac mass m is only one of many equivalent ways of breaking 
this symmetry down to a t/(l) x t/(l) subgroup. Indeed, by applying a general 
U (2) transformation, we find that the Dirac mass term can be turned into a linear 
combination of the following three mass terms: m, y 3 ?^, and iy 5 ms. (In principle, 
there is also a possibility of the so-called Haldane mass term y i "f'A, which is a 
singlet under U(2). We do not discuss it here. However, as we will see in Sec. 13.31 
the Haldane mass plays an important role in graphene.) 

In our perturbative analysis, we did not get any nonzero or m$ because the 
vacuum alignment was predetermined by a "seed" Dirac mass mo in the free theory, 
see Eq. dS). 



2.2.2 Large N approximation 

It is instructive to generalize the above analysis in the NJL model to the case of 
strong coupling. While magnetic catalysis occurs even at arbitrarily weak coupling, 
such a generalization will be useful to understand how magnetic catalysis is lost in 
the limit of the vanishing magnetic field. 
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At strong coupling, a reliable solution to the NJL model can be obtained by using 
the so-called large N approximation, which is rigorously justified when the fermion 
fields in Eq. (0 carry an additional, "color" index a = l,2,...,N, and N is large. 
Using the Hubbard-Stratonovich transformation 1178111151 . one can show that the 
NJL theory in Eq. (O is equivalent to the following one: 



if = Ti>f D„W - <¥{o + fx + ifn)W- 



1 

2G 



(a 2 + 7t 2 + T 2 ) 



Note that the equations of motion for the new composite fields read 
a = -G(WW), t= -G^f^), n = -G{Wif x ¥). 



(13) 



(14) 



Under U (2) flavor symmetry transformations, these composite fields transform into 
linear combinations of one another, but the quantity a 2 + % 2 + x 2 remains invariant. 



The effective action for the composite fields, 



1 

2G 



d 3 x(a 2 ■ 



% ) — /trLn 



iy^D/i - (a + fx + ifn) 



(15) 



is obtained by integrating out the fermionic degrees of freedom from the action. It is 
convenient to expand this effective action in powers of derivatives of the composite 
fields. The leading order in such an expansion is the effective potential V (up to the 
minus sign). Because of the flavor symmetry, the effective potential depends on a, 
%, and x fields only through their 1/(2) invariant combination p 2 = a 2 + n 2 + x 2 . 

Using the proper-time regularization, one obtains the following explicit expres- 
sion for the effective potential l99l 11021 : 



v ( p)4 

n 



2\eB\ 



p\eB\ 



, (16) 



where we dropped the terms suppressed by the ultraviolet cutoff parameter A and 
introduced a dimensionless coupling constant g = NAG/n. 

The field configuration p that minimizes the effective potential is determined by 
solving the equation dV /dp = 0, i.e., 



.1 1 \ \eB\ \eB\ ( 1 , p l 



(17) 



In essence, this is the gap equation. At weak coupling, g — > 0, in particular, we 
obtain the following approximate solution: 



m = p n 



GN- 



\eB\ 
2k 



(18) 



which is the large N generalization of the result for the dynamical mass in Eq. (fT2b . 
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2.2.3 Zero magnetic field limit in 2 + 1 dimensions 

Before concluding the discussion of the (2 + 1) -dimensional NJL model, it is in- 
structive to consider how the above analysis of flavor symmetry breaking modifies 
in the zero magnetic field case. By taking the limit B — >• in Eq. (fTol i. we arrive at 
the following effective potential: 

N 

y B=0 (p) * - 

K 

Now, the zero-field limit of the gap equation, dVB=o/dp = 0, is given by 

A(-L-I)p=p=. (20) 

This is very different from Eq. ( TPTl i. In particular, the only solution to this equation 
at g < y/li is a trivial one, p = 0. The nontrivial solution appears only in the case of 
a sufficiently strong coupling constant, g > y/n. This is in a stark contrast with the 
dynamical mass generation in the presence of a magnetic field, where a nontrivial 
solution exists at arbitrarily small values of the coupling constant g. 



A /l 1 \ 2 1 3 



(19) 



2.3 Magnetic catalysis in 3+1 dimensions 

Let us now extend the analysis of the previous subsection to the case of a (3 + 1)- 
dimensional model, where the dynamics is truly nonperturbative. The Lagrangian 
density of the corresponding NJL model reads 



(21) 



This model possesses the U(1)l x U(1)r chiral symmetry. The symmetry will be 
spontaneously broken down to the U ( 1 )l+r subgroup when a dynamical Dirac mass 
is generated. 



2.3.1 Weak coupling approximation 

In the weakly coupled limit, the gap equation in the mean-field approximation reads 

m = Gtr[S(x,x)]. (22) 

Formally, it is same as the gap equation in the (2 + 1 ) -dimensional model in Eq. (fTTT ). 
However, as we will see now, its symmetry breaking solution will be qualitatively 
different. 



Magnetic Catalysis: A Review 



11 



Let us start by showing that the chiral condensate vanishes in the free theory in 
3 + 1 dimensions when the bare mass goes to zero, mo 0. By making use of the 
proper-time representation, see Eq. (|97| i in the Appendix, we obtain 



(WW) 



-tr[5(x,x)] 



m \eB\ 



mp 
(2k) 2 



A" —mi In 



[2n) 2 
A 2 



A- 



leSlln- 



2\eB 
\eB\ 



Ttm, 



ds 
— t 

l/A 2 S 

-Ye 
- mo In 



"o coth(i|eB| 
,2 



I d?| to- 



rn 



\%\eB\ 



2|eB|lnr 



2\eB\ 



A- 



o 



2eB 



O 



\eB\ 



(23) 



In the limit mo — > 0, we see that the condensate indeed vanishes. This means that 
we cannot apply a perturbative approach to find any nontrivial (symmetry breaking) 
solutions to the gap equation in Eq. ( l22l . 

The explicit form of the gap equation ( 1221 reads 



{2nf 

Its nontrivial solution is given by 



eB\ 



■ exp 



A" 



A^ 



leBlln- 



eB 



7Tm z 



2eB 



exp 



2n 2 
G\eB\ 



(24) 



(25) 



When G — > this result reveals an essential singularity. Obviously, such a depen- 
dence cannot possibly be obtained by resuming any finite number of perturbative 
corrections in powers of a small coupling constant. Therefore, despite the weak 
coupling, the result for the dynamical mass (l25l l is truly nonperturbative. 



2.3.2 Zero magnetic field limit in 3 + 1 dimensions 

It is instructive to compare the above dynamics of spontaneous symmetry breaking 
with case of the zero magnetic field. At B = 0, the chiral condensate in the free 
theory is easily obtained by taking the appropriate limit in Eq. ( 1231 . i.e., 



(WW) 



/)?o 



B=0 



{2Kf 

The corresponding gap equation is 



A- 



ml I ln^r + 1 -y E 



(26) 



(2K) 



A- 



2, A 

m \n^r 



(27) 



Because of the negative sign in front of the logarithmic term, this equation does not 
have any nontrivial solutions for the dynamical mass at vanishingly small coupling 
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constant g = GA 2 /(2n) 2 . In fact, for the whole range of subcritical values, g < 1, 
the only solution to this gap equation is m = 0. The nontrivial solution appears only 
in the case of sufficiently strong coupling, g > 1. 



2.4 Symmetry breaking as bound state problem 

In this subsection, we consider an alternative approach to the problem of chiral 
symmetry breaking in the NJL model in a constant magnetic field. As we will see, 
this approach is particularly beneficial for illuminating the role of the dimensional 
reduction in magnetic catalysis. 

Instead of solving the gap equation, we consider the problem of bound states with 
the quantum numbers of the Nambu-Goldstone bosons, using the method of a homo- 
geneous Bethe-Salpeter equation (for a review, see Ref. IU42ID . The underlying idea 
for this framework is motivated by the Goldstone theorem [148, 82, 83 1. The theo- 
rem states that spontaneous breaking of a continuous global symmetry leads to the 
appearance of new massless scalar particles (i.e., Nambu-Goldstone bosons) in the 
low-energy spectrum of the theory. The total number of Nambu-Goldstone bosons 
and their quantum numbers are determined by the broken symmetry generators. 

The homogeneous Bethe-Salpeter equation for a pion-like state takes the form 
0142111031 : 

Xab(x,T,P) = -i J d A x'd A y'd A x"d A y"S im (xj)^^ (*'/,*"/') 

x Xa 2 b 2 (x"y;P)S blh {y',y), (28) 

where Xab{x,y,P) = (0\T\l/ a (x)\j/b(y)\P;it) is the Bethe-Salpeter wave function of 
the bound state boson with four-momentum P, and S c ,h(x,y) = (Q\T\j/ a (x)\jfb(y)\0) 
is the fermion propagator. Here and below, the sum over repeated Dirac indices (a\, 
b\, «2, bj) is assumed. The explicit form the Bethe-Salpeter kernel is II142II1031 : 

Ka l b l ;a 2 b 2 (x'y',x"y") = G [Kbi 8 b 2 a 2 + (*1§)*2«a2 ~ ('ft)aia 2 (^)*2*1 

- Ken Kb, ] S 4 (x' -y')8 4 (x' - x" ) 8 4 (x' -y"). (29) 

It is convenient to rewrite the wave function in terms of the relative coordinate z = 
x — y and the center of mass coordinate X = (x+y)/2, 

Xa b (X,z;P) = e is ^/ l2 e- iP » x *% ab (X,z;P), (30) 

where we factorized the Schwinger phase factor, see Eq. d88l . and introduced the 
notation: s± = sign(efi) and l=\j \J\eB\. After substituting the wave function < f30b 
and the kernel d29b into Eq. d28l i, we arrive at the following equation: 



Xab(z;P) = -iG J d 4 X'S aai (| -X 1 ) 



Kh tr[z(0;P)]-(7 5 W 1 tr[r 5 z(0; J P)] 
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-Xa 1 b 1 (0;P) + (j 5 ) ai a 2 Xa 2 b 2 {0;P)(75)b 2 b l 



Sb,b ( 2 + x ' 



x ^(^)/^. (31) 

Here we took into account that the equation admits a translation invariant solution 
and replaced Xab(X,z;P) — > Xab(z,P). Note that, on the right-hand side of Eq. (|3~TT ), 
the dependence on the center of mass coordinate X completely disappeared after a 
shift of the integration variable X' — > X' — X was made. 

In the lowest Landau level approximation, one can show that the Fourier trans- 
form of the Bethe-Salpeter wave function takes the form 01031 : 

JiaCp^-X) )=A(p|| e "J / ' /' J, (32) 

where = (fi),p 3 ),p^ = (p 1 ) 2 + {p 2 ) 2 , and ^± = \ (l ±/sxy'y 2 )- The new func- 
tion A(/»||) satisfies the equation: 

. G|efl| [ A(k\\)d\ 

A(/, " )= ^y^w' (33) 

where we made the Wick rotation (o — > j'eo). The solution to this equation is a 
constant: A(pii) = C. Dropping nonzero C and cutting off the integration at A, we 
finally arrive at the gap equation for the mass parameter m: 

G\eB\ [ Al dk 2 

(34) 



4-n 2 Jo k 2 + m 2 
The solution to this equation is 



w " A2exp lrcpiJ' (35) 

which, to leading order, agrees with the solution obtained in Eq. ( f25b . 

The Bethe-Salpeter equation ( |33l can be rewritten in the form of a two-dimen- 
sional Schrodinger equation with an attractive 5-function potential. In order to see 
this explicitly, let us introduce the following wave function 



2% 2 



d 2 k\\ --'' k " r 
(2nj 2 k 2 +m 2 



¥(r)= I fn J 2a ,_, A(k\\). (36) 



Taking into account that A(p») satisfies Eq. d33l ). it is straightforward to show that 
the wave function \jf(r) satisfies the following Schrodinger type equation: 



+ m 2_^i 5 2 (r))v , (r) = 0) 
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in which — m 2 plays the role of the energy E. Since m 2 must be positive, the prob- 
lem is reduced to finding the spectrum of bound states (with E = —m 2 < 0) in the 
Schrodinger problem. The potential energy in Eq. ( f37T > is expressed in terms of 

~ r d 2 k\\ , , 

«< r K WF' ' ' <38) 

which is a regularized version of the 5-function that describes the local interaction 
in the NJL model. 

Notice that, by using the same approach, one can show that the Bethe-Salpeter 
equation for a massless NG-boson state in the NJL model in 2 + 1 dimensions can 
be reduced to the gap equation 

Mp ) -9M[ A iW* (39) 
A{P) ~ 2k 2 J- A k 2 +m 2 - (39) 

This has the same solution for the mass as in Eq. (fT2l . Also, this integral equation 
is equivalent to the following one-dimensional Schrodinger equation: 

d 2 , G\eB\ „ . A , s 
-^2+m 2 - -^$A to J Vto = 0, (40) 

where the regularized version of the 5-function is given by 5 A (x) = J_ A j^e~ lkx . 



2.5 Analogy with superconductivity 

It is interesting to point that the dynamics described by the gap equation in the case 
of magnetic catalysis has a lot of conceptual similarities to the mechanism of super- 
conductivity in metals and alloys. This is despite the clear differences between the 
two phenomena that we discussed in the Introduction. (In order to avoid a possible 
confusion, let us emphasize that here we compare the nonrelativistic Cooper pairing 
dynamics in superconductivity in the absence of magnetic fields with the relativistic 
dynamical generation of a mass in the presence of a constant magnetic field.) 

The corresponding gap equation in the Bardeen-Cooper-Schrieffer theory ifTTl of 
superconductivity can be written in the following form: 

rtia>D de 

l=GN(0) _ __— , (41) 

Jo Ve 2 + A 2 

where N(Q) is the density of electron states at the Fermi surface, (Od is the Debye 
frequency, and A is the energy gap associated with superconductivity. The solution 
for the gap reads 

1 



4 ~ftflto«p [-__). (42) 
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At weak coupling, this solution has the same essential singularity as the dynamical 
mass parameter in Eq. d35l l. We can argue that the similarity is not accidental. To 
see this clearly, let us rewrite the gap equation in the problem of magnetic catalysis 
in the lowest Landau level approximation, see Eq. ( f34t . as follows: 



where the Wick rotation was performed (to — > ico). As we see, the structure of 
this gap equation is identical to its counterpart in the BCS theory after we iden- 
tify the density of states N(0) with the density of states in the lowest Landau level, 
\eB\/ (2k), and the Debye frequency (Od with the cutoff parameter A . 

The similarity between the BCS theory of superconductivity and magnetic catal- 
ysis goes deeper. In particular, the generation of a nonzero gap in superconductors 
can be also thought of as the result ofa3 + l^l + l dimensional reduction of 
the phase space around the Fermi surface. Also, just like in magnetic catalysis, it is 
essential that the density of states at the Fermi surface is nonzero. 



2.6 Bound states in lower dimensions 

As we saw in Sec. 12.41 the problem of spontaneous symmetry breaking and the asso- 
ciated dynamical generation of the Dirac mass can be reformulated as a problem of 
composite massless states with the quantum numbers of Nambu-Goldstone bosons. 

In the presence of a constant magnetic field, in particular, we also found that 
the corresponding Bethe-Salpeter equation for the bound states can be recast in an 
equivalent form as a Schrodinger equation in a dimensionally reduced space. The 
dimensional reduction is D — > D — 2 and, therefore, the relevant problem of bound 
states is considered in spaces of lower dimensions. 

In order to prove that the essence of magnetic catalysis is directly connected with 
this reduction, let us consider a simple quantum mechanical problem: the formation 
of bound states in a shallow potential well in spaces of various dimensions. As 
we will see, at least one bound state does exist in one- and two-dimensional cases 
M175l|176l[l9lll34| , but not always in three dimensions. We will also see that, while 
the result for the binding energy is perturbative in the coupling constant in one 
dimension, it has an essential singularity in two dimensions. 

2.6.1 Bound states in a one-dimensional potential well 

Let us start from the simplest one-dimensional problem of a nonrelativistic particle 
of mass to* confined to move on a line. Let the potential energy of the well be given 
by U(x), which is negative and quickly approaches zero when |jc| — s- °°. One can 
show that even a vanishingly small depth of the potential well is sufficient to pro- 




(43) 
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duce a bound state (i.e., a quantum state with a negative energy). The corresponding 
binding energy is given by II 1341 



If we rescale the potential energy U(x) by a "coupling constant" factor g, i.e., 
U (x) — > gU(x), we find that |£id| ~ g 2 as g — > 0. In other words, the binding energy 
has a power-law dependence as a function of the depth of the potential energy U (x). 
This is a typical result that can be obtained by perturbative techniques, controlled 
by powers of the small parameter g |fl9l . 

The above conclusion remains valid basically for any attractive potential U(x). 
For example, one can rigorously prove that, if J (I + \x\)\U(x)\dx < °°, there is a 
bound state for all small positive g if and only if J U(x)dx < (i.e., the potential is 
attractive at least on average) 111251 . 

2.6.2 Bound states in a two-dimensional potential well 

In the case of a two-dimensional system (i.e., a nonrelativistic particle of mass 
confined to move on a plane), the general conclusion about the existence of a bound 
state around a potential well of a vanishingly small depth still remains valid. How- 
ever, an important qualitative difference appears in the result. The binding energy 
reveals an essential singularity as a function of the depth of the potential well. In or- 
der to understand this better, let us consider a problem with a cylindrically symmet- 
ric potential energy U(r), where r is the radial polar coordinate in the plane. If the 
potential energy is sufficiently shallow and localized (i.e., | rU(r)dr\ <C m t /h 2 ), 
one finds that the energy of the bound state is given by II 1341 



where a is the characteristic size of the potential well. The fact that this energy is 
singular can be made explicit by rescaling the potential energy U (r): U(r) — > gU (r). 
Then, we find that |2?2d| ~ e x P( — C/g) as g — > (here C is a constant determined by 
the shape of the potential well). Unlike the g 2 power-law suppression of the bind- 
ing energy in one dimension, this is a much stronger suppression indicating a much 
weaker binding. Moreover and perhaps more importantly, such an essential singu- 
larity cannot possibly be obtained by resuming any finite number of perturbative 
corrections, controlled by powers of the small parameter g. Therefore, the singular 
behavior of the binding energy in two dimensions is a sign of a truly nonperturbative 
(albeit weakly-interacting) physics. 

Again, this result is very general. It can be rigorously proven that, in the case 
when / \U(x)\ l+£ d 2 x < °o (with some e > 0) and J (I +x 2 ) £ \U(x)\d 2 x < °°, there is 




(44) 




(45) 
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a bound state for all small positive g if and only if / U (x)d 2 x < (i.e., the potential 
is attractive at least on average) II175II1761 . 



2.6.3 Bound states in a three-dimensional potential well 

Now, in the three-dimensional case, there are no bound states if the potential well 
is too shallow in depth. This was first shown by Peierls in 1929 11581 . This can be 
demonstrated, for example, in a special case of a spherically symmetric potential 
well of a finite size, 

r for r<a, 



U(r) = i 6 8^? - ' (46) 

for r > a. 



The condition to have at least one bound state is g > 1 113411 . In other words, the 
depths of the potential well (or the strength of the "coupling constant" g) should be 
larger than the critical value, given by g cr = 1 . In the supercritical regime, g = 1 + e 
with < e <C 1, the binding energy is given by 01341 

Nl = ^ 2 (47) 
In the subcritical regime g < 1, on the other hand, there are no bound states at all. 



3 Magnetic Catalysis in Gauge Theories 



Motivated by the fact that magnetic catalysis has a rather general underlying 
physics, explained by the dimensional reduction of the particle-antiparticle pairing, 
it is natural to ask how it is realized in gauge theories with long-range interactions, 
such as QED. This problem was discussed in numerous studies IU01lll36lll56|[T57l 
[IIO][I(I]|5l][l35][I05]|^ Here we will briefly 

review only the key results and refer the reader to the original papers for further 
details. 



3.1 Magnetic Catalysis in QED 

Using the same conceptual approach as outlined in Sec. 12.41 for the NJL model, 
one can show that, in Euclidean space, the equation describing a pion-like Nambu- 
Goldstone boson in QED in a magnetic field has the form of a two-dimensional 
Schrodinger equation 11011 : 
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^2- — +™ 2 + V(r) 



d- 2 



V/(r)=0. 



(48) 



The function yf(r) is defined in terms of the Bethe-Salpeter wave function A(p) in 
exactly the same way as in the NJL model, see Eq. (l36l l. This time, however, A(p) 
satisfies a different integral equation, 



A(p) 



a 

2K 2 



d 2 kA(k) f dx exp(-x£ 2 /2) 



k 2 



(k-p) 



(49) 



where l=\j \J\eB\ is the magnetic length. Note that, in addition to using the low- 
est Landau level approximation, we assumed that the photon screening effects are 
negligible. As is easy to check, the explicit form of the potential V(r) is given by 

2l 2 



V{r) = W eXP {2f 2 



(50) 



where r 2 = r\ + ri and Ei(x) = — f™ x dtexp(— t)/t is the integral exponential func- 
tion l94l . Since V (r) is negative, we have a Schrodinger equation with an attractive 
potential, in which the parameter — m 2 plays the role of the energy E. Therefore, the 
problem is again reduced to finding the spectrum of bound states with E = —m 2 < 0. 

It is known that the energy of the lowest level E(a) for the two-dimensional 
Schrodinger equation is a nonanalytic function of the coupling constant a at a = 
111761 . If the potential V(r) were short-range, the result would have the form m 2 = 
—E(a) <x exp[— l/(Ca)], where C is a positive constant 11751 11761 . In our case, 
however, we have a long-range potential. Indeed, using the asymptotic expansion 
forEi(x) lH, we get: 



2a 1 



7t r z 



(51) 



In order to find an approximate solution for m 2 , one can use the integral equation 
(l49b at p — 0. As a — >• 0, the dominant contribution in the integral on the right-hand 
side comes from the infrared region k 2 < m 2 . Therefore, 



a , . f d 2 k f dxexp(— y/2) a 
A(0 ) ~ A (0) / I J} "! ' * -A(0) 



k 2 + m 2 







l 2 k 2 +y 



4-n 



In 



V- 



(52) 



which implies that 11011 



(53) 



A slightly more careful analysis of the integral equation d49b can be made by ap- 
proximating the interaction kernel so that the exchange momentum (k — p) 2 in the 
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denominator is replaced by ma\(k 2 ,p 2 ). The problem then reduces to an ordinary 
differential equation with two (infrared and ultraviolet) boundary conditions. The 
approximate analytical solution reveals that the lowest energy bound state, which 
describes the stable vacuum solution in quantum field theory, corresponds to the 
following value of the dynamical mass 11011 : 



Unfortunately, the approximation used in this analysis is not completely reliable. 
There are higher order diagrams that can substantially modify the interaction poten- 
tial and, in turn, the result for the dynamical mass. For example, taking into account 
the vacuum polarization effects in the improved rainbow (ladder) approximation, in 
which the free photon propagator is replaced by a screened interaction with the one- 
loop photon self-energy, the result changes. The corrected expression for the mass 
has the same form as in Eq. ( f54b . but with a replaced by a/2 11011 . This is a clear 
indication that, despite weak coupling, there can exist other relevant contributions, 
coming from higher order diagrams. 

A further study showed that, by using a similarity between the magnetic catalysis 
problem in QED and the exactly solvable Schwinger model H167ll73l . one can find a 
special nonlocal gauge, in which the leading singularity of the dynamical mass can 
be extracted exactly 11041 , 



where N is the number of fermion flavors, F(a) ~ (Na) 1 / 3 , C\ ~ 1.82 ±0.06 and 
C ~ 0(1). Note that the leading singularity in the final expression for the mass is 
quite different from that in the rainbow approximation ( |54| ). 

The magnetic catalysis of chiral symmetry breaking in QED yields a rare exam- 
ple of dynamical symmetry breaking in a (3 + 1) -dimensional gauge theory with- 
out fundamental scalar fields, in which there exists a consistent truncation of the 
Schwinger-Dyson equation. 



3.2 Magnetic catalysis in QCD 

Recently there was an increased interest in studies of QCD in a strong magnetic field 



imiiiMiiiiiiainaiEiissi^^ 



There are several reasons why such investigations may be of interest. Very strong 
magnetic fields are known to have existed in the Early Universe |181|l44l[30l[T8ll95l 
and are expected to be generated in relativistic heavy ion collisions 1122111771 . Since 
the chiral symmetry plays a profound role in QCD, it is interesting to study also the 
role of magnetic catalysis in this theory 11441112111 . 




(54) 




(55) 
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Because of the property of asymptotic freedom, one can argue that the dynamics 
underlying magnetic catalysis in QCD is, at least in principle, weakly coupled at 
sufficiently large magnetic fields 11211 . This fact can be used to justify a consistent 
truncation of the Schwinger-Dyson equation, resembling that in QED, which we 
discussed in the preceding section. 

Let us start by introducing a QCD like theory with N u up flavors of quarks having 
electric charges 2e/3 and N d down flavors of quarks having electric charges —e/3. 
(The total number of flavors is Nf = N u +N d .) It is important to distinguish the up 
and down types of quarks because the chiral symmetry subgroup that mixes them 
is explicitly broken by the external magnetic field. Taking this into account, we find 
that the model is invariant under the SU(N U ) L x SU(N U ) R x SU(N d ) L x SU(N d ) R x 
U^\\)a chiral symmetry. The anomaly free subgroup U^'(1)a is connected with 
the conserved current which is the difference of the U^\1)a and 11^(1) a currents. 
[The U(>(1)a symmetry is of course absent when either N d or N u equals zero.] 
A dynamical generation of quark masses spontaneously breaks the chiral symme- 
try down to SU {Nu)v x SU(N ( /)v and gives rise to N% + N% — 1 massless Nambu- 
Goldstone bosons in the low-energy spectrum. 

Just like in QED, the vacuum polarization effects play a very important role in 
QCD in the presence of a strong magnetic field. By properly modifying the known 
result from the Abelian gauge theory IU38l[36ll26l to the case of QCD, we find that 
the gluon polarization tensor has the following behavior: 

^ v ^ 5AB te-« v )l^ for i*ii«^> ^ 

n AB.,v „ _^ 8 AB ^ k v _ k 2 g ^ £ Ml i for m 2 <<; ^ <<; ^ (5?) 

where = g^kv and g^ v = diag( 1 , 0, 0, — 1 ) is the projector onto the longitudinal 
subspace. Notice that quarks in a strong magnetic field do not couple to the trans- 
verse subspace spanned by g^ = g^ v — g^ v = diag(0, — 1, — 1,0) and k^ = g^ky. 
This is connected with the dominant role of the lowest Landau level, in which quarks 
are polarized along the magnetic field. 

The expressions d56l l and d57l i coincide with those for the polarization operator 
in the massive Schwinger model II 1671 if the parameter a s \e q B\/2 here is replaced 
by the dimensional coupling a,\ of (1 + l)-dimensional QED. In particular, Eq. d57l i 
implies that there is a massive gluon resonance with the mass given by 

N f 

(if a 

M l = E -Ml = (2N u +N d )^\eB\. (58) 

q=\ U iK 

This is reminiscent of the pseudo-Higgs effect in the (l+l)-dimensional massive 
QED. It is not the genuine Higgs effect because there is no complete screening of 
gluons in the far infrared region with |&|| <§; m^, see Eq. (156) . Nevertheless, the 
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pseudo-Higgs effect is manifested in creating a massive resonance and this reso- 
nance provides the dominant force leading to chiral symmetry breaking. 

In the end, the dynamics in QCD in a strong magnetic field appears to be essen- 
tially the same as in QED, except for purely kinematic changes. After expressing 
the magnetic field in terms of the running coupling a s at the scale \/|c5| using 



l^ln leB| 



a. 



A 



where b = 



QCD 



1 IN C - 2N f 
12n ! 



we obtain the result for the dynamical mass in the following form 11441 : 



2Ci 



A2CD(c 9 a 4 .) 2/3 exp 



ba s a s (N? - \)\n(C 2 /c q a s ) 



(59) 



(60) 



where e q is the electric charge of the q-th quark and N c is the number of colors. The 
numerical factors C\ and C2 are of order 1, and the value of c q is given by 



1 

67T 



(2N U +N d ) 



(61) 



Because of the difference in electric charges, the dynamical mass of the up-type 
quarks is considerably larger than that of the down-type quarks. 

It is interesting to point that the dynamical quark masses in a wide range of 
strong magnetic fields, Aq CD <C \eB\ < (10 TeV) 2 , remain much smaller than the 

dynamical (constituent) masses of quarks mf ] -300 MeV in vacuum QCD without 
a magnetic field. This may suggest that QCD can have an intermediate regime, in 
which the magnetic field is strong enough to provide a gluon screening to interfere 
with the vacuum pairing dynamics [144 76 1, but not sufficiently strong to produce 
large dynamical masses through magnetic catalysis. In this intermediate regime, the 
dynamical mass and the associated chiral condensate could be decreasing with the 
magnetic field. The corresponding regime may start already at magnetic fields as 
low as 10 I9 G, when the gluon mass M g , given by Eq. (l58l l. becomes comparable to 
Aqcd- (For me estimate, we assumed that the value of the coupling constant is of 
order 1 at the QCD energy scale.) 



3.3 Magnetic Catalysis in Graphene 

In this section, we briefly discuss the application of the magnetic catalysis ideas to 
graphene in the regime of the quantum Hall effect. 

Graphene is a single atomic layer of graphite 11541 that has many interesting 
properties and promises widespread applications (for reviews, see Refs. HI 091 l27l 
[l]). The uniqueness of graphene is largely due to its unusual band structure with 
two Dirac points at the corners of the Brillouin zone. Its low-energy excitations are 
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described by massless Dirac fermions H170II . Because of a relatively small Fermi ve- 
locity of quasiparticles, vp ~ c/300, the effecting coupling constant for the Coulomb 
interaction in graphene, a = e 2 /( e v F )■> is about 300 times larger than the fine struc- 
ture constant in QED, e 2 / (fioc) ~ 1 /137. 

When graphene is placed in a perpendicular magnetic field, it reveals an anoma- 
lous quantum Hall effect |[T53l[TF7l , exactly as predicted in theory Ifl89l [1081 [1591 . 
The anomalous plateaus in the Hall conductivity are observed at the filling factors 
v = ±4(« + 1 /2), where n = 0, 1,2, . . . is the Landau level index. The factor 4 in 
the filling factor is due to a fourfold (spin and valley) degeneracy of each Landau 
level. As for the half-integer shift in the filling factor, it is directly connected with 
the Dirac nature of quasiparticles ifTTOl \U2\ [1241 [1231 l88l . 

It was observed experimentally Ml 861 11201 [37l l22l that there appear additional 
plateaus in the Hall conductivity when graphene is placed in a very strong mag- 
netic field. The new plateaus can be interpreted as the result of lifting the four- 
fold degeneracy of the Landau levels. In the case of the lowest Landau level, in 
particular, some of the degeneracy, i.e., between the particle and hole states, can 
be removed when there is a dynamical generation of a Dirac mass. Considering 
the possibility of magnetic catalysis, such an outcome seems almost unavoidable 



The low-energy quasiparticle excitations in graphene are conveniently described 
in terms of four-component Dirac spinors = {^kas,Wkbs,Wk'BsiWk'As)^ intro- 
duced for each spin state s =t,|- Note that the components of f v are the Bloch 
states from two sublattices (A,B) of the graphene hexagonal lattice and two valleys 
(K, K') at the opposite corners of the Brillouin zone. The approximate low-energy 
Hamiltonian, including the kinetic and Coulomb interaction terms, is given by 



where Uc(r) is the Coulomb potential, which takes into account the polarization 
effects in a magnetic field l88l l9ll . Note that the two electron spins (s ~f,i) in 
graphene give rise to two independent species of Dirac fermions. As a result, the 
Hamiltonian possesses an approximate U(4) symmetry 11881 . which is a generaliza- 
tion of the U(2) flavor symmetry discussed in the case of the one-species model in 
Sec. 12.21 The 16 generators of the extended U(4) flavor symmetry are obtained by 
a direct product of the 4 generators of the U(2) group acting in the valley space 
(K,K'), and the 4 generators of the U (2) spin symmetry. 

The U (4) symmetry is preserved even when the electron chemical potential term, 
— /if ^f, is added. The inclusion of the Zeeman term, which distinguishes the elec- 
tron states with opposite spins, breaks the symmetry down to the U^(2) x t/|(2) 
subgroup. The explicit form of the Zeeman term is given by /igBf ^cr^f , where B 
is the magnetic field, = eh/{2mc) is the Bohr magneton, and (73 is the third Pauli 
matrix in spin space. An interesting thing is that this explicit symmetry breaking is 



H22[T23][88 98 89l fTT3l fl72l I9l1 . 




(62) 
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a small effect even in very strong magnetic fields. To see this, we can compare the 
Zeeman energy Ez with the Landau energy Eg , 

£ Z = ld B B = 5.8 x l(T 2 fl[T] meV, (63) 
Ei = ^Tiv 2 F \eB\/c = 26 y/~B\J] meV. (64) 

Therefore, the Zeeman energy is less then a few percent of the Landau energy even 
for the largest (continuous) magnetic fields created in a laboratory, B < 50 T. 

Because of the large flavor symmetry, there are many potential ways how it can 
be broken l89l I9TI . Here we mention only the possibilities that are connected to 
the magnetic catalysis scenario at zero filling v = (i.e., the lowest Landau level is 
half-filled). 

We will allow independent symmetry breaking condensates for fermions with 
opposite spins. Also, in addition to the usual condensates (no sum over the 

repeated spin indices here), we introduce the time reversal odd ones, 
ll89ll9T1l . While the former will give rise to Dirac masses m s (s =t,4) in the low- 
energy theory, the latter will result in the Haldane masses A s {s =t,4) 01 121 . 

In the ground state, one can also have additional condensates, (f ^<7 3 f ) and 
(fty^P^f), associated with nonzero spin and pseudo-spin (valley) densities. To 
capture this possibility in the variational ansatz, one needs to include a spin chemical 
potential /I3 and two pseudo-spin chemical potentials p, s (s =t,4)- Thus, the general 
structure of the (inverse) full fermion propagator for quasiparticles of a fixed spin 
has the following form: 

S;\co;r,r') = -i[fco-v F (JC-y) + E+] 8 2 (r-r'), (65) 

where the generalized self-energy operator E + is given by 

E + = -m s + yVv + isj-fffis + is±fy L y 1 A s . (66) 

Functions m s , fi s , p, s , and A s on the right-hand side depend on the operator valued 
argument [it ■ y) 2 £ 2 , whose eigenvalues are nonpositive even integers: — 2m, where 
n = 0, 1,2, . . .. Therefore, in the Landau level representation, m s , jJ. s , fi. s , and A s will 
get an additional Landau index n dependence: m n i , fi n s , p, n s , and A„ , s . 

The Schwinger-Dyson equation for the full fermion propagator takes the form 

r 1 (r-(';r,r')=X l ((-(';r,r')+eYx('-'';r, r')fD{t' -t;r'-r), (67) 

where D(t;r) is the photon propagator mediating the Coulomb interaction. The lat- 
ter is approximately instantaneous because the quasiparticle velocities are much 
smaller than the speed of light. In momentum space, the photon propagator takes 
the following form: 
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where n(0,k) is the static polarization function and £o is a dielectric constant. 

It should be noted that, in the coordinate-space representation, both the fermion 
propagator and its inverse contain exactly the same Schwinger phase, see Eq. (|88l ). 
After omitting such a (nonzero) phase on both sides of Eq. d67l i and performing the 
Fourier transform with respect to the time variable, we will arrive at the following 
equation for the translationally invariant part of the fermion propagator l9D : 

rw>=V(-) + |/:^f§^y>^>A m 

In the Landau level representation, this equation is equivalent to a coupled set of 
4 x 2 x « max equations, where we counted 4 parameters (m, jj,, fi, and A), 2 spins 
(s =f, I), and n max ~ [A 2 / (2 \eB | )] Landau levels below the ultraviolet energy cutoff 
A, where the low-energy theory is valid. 

The explicit form of the gap equations can be found elsewhere |9"T1 . The corre- 
sponding set of equations can be solved by making use of numerical methods. Here, 
instead, we will discuss only some general features of the solutions in the lowest 
Landau level approximation, which can be obtained with analytical methods. 

Let us start by considering the solutions to the gap equations for quasiparticles 
of a fixed spin. In the lowest Landau level approximation, there are two independent 
gap equations, i.e., 

Meff - M = -J- ^0 [»F (m e ff ~ Meff ) - «F Oeff + Meff)] i (70) 

m e s = [1 - n F Oeff - Meff) - n F (m eS + jLi e ff)] , (71) 

where a = e 2 / (£oVf) rj 2.2/ Eq is the coupling constant, np{x) = 1/ (e*/ T + l) is 
the Fermi distribution function, and J£q lS the interaction kernel due to the Coulomb 
interaction in the lowest Landau level approximation. In the above equations, we 
used the shorthand notation jLi e ff = /J. — A and m e ff = m — p. for the two independent 
combination of parameters that determine the spectrum of the lowest Landau level 
quasiparticles, 

CO- = -jLieff - meg, and co + = -^ eff + m ef[ . (72) 
At zero temperature, the gap equations reduce down to 

CC£f 

Meff = M + — 7K=sign(jUeff)0(|Aieff I - Kff|), (73) 
4V27T 

\m e g\ = -^L(?(|m eff | - | Meff |). (74) 
4V27T 

Here we used the value for the interaction kernel J£q = 1 / {2\f2n), which is obtained 
in the approximation with screening effects neglected 1911 . One of the solutions to 
this set of equations has a nonzero dynamical Dirac mass (m °= OCQ), i.e., 
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«££ n ae e ae e 
m eff =_^, 4=0, jl=<li<—F=- (75) 



The other two solutions have nonzero Haldane masses (A <*; CCEg), i.e., 

m eff = 0, 4 = —^, _oo< M <_' (76) 
4V27T 4V2^ 



m eff = 0, 4 = -—^=, -— , =L<Ai<-. (77) 

In both types of solutions, the values of the masses are proportional to a power of the 
coupling constant a, as expected from the dimensional reduction 1113411175111761 . 

In order to determine the ground state in graphene when both spin states are 
accounted for, one has to find among many possible solutions the one with the lowest 
free energy. In the approximation used here, the ground state solution at v = filling 
(i.e., an analog of the vacuum state in particle physics) corresponds to a spin-singlet 
state with equal in magnitude, but opposite in sign Haldane masses for the two spin 
states [|9 11 : Af = —A^, i.e., a mixture of the two solutions in Eqs. (l76*l l and ( [77) . 

The symmetry of the corresponding ground state is t/f (2) x t/j.(2), but with the 
Zeeman energy splitting dynamically enhanced by the nonzero Haldane masses. The 
quasiparticle energies of the dynamically modified lowest Landau level are |9"T1 

co r = -n + £ Z +\A t \>0, (x2), (78) 
© ; = -M-£z-|4;|<0, (x2), (79) 

which show that the original fourfold degeneracy is indeed partially lifted. 



4 Concluding Remarks 

We hope that this review of magnetic catalysis is sufficient to convey the main idea 
of the phenomenon in terms of simple and rather general physics concepts. From 
the outset, this review was never meant to be comprehensive. Here we concentrated 
only on the bare minimum needed to understand the phenomenon as a consequence 
of the underlying dimensional reduction of the fermion-antifermion pairing in a 
magnetic field [99 , 1 1 021 1 1 001 1 1 ll . For further reading and for deeper insights into 
various aspects of the magnetic catalysis, it is suggested that the reader refers to the 
original literature on the topic. 

Over nearly 20 years of research, there has been a lot of progress made in our 
understanding of magnetic catalysis. A rather long list of research papers at the end 
of this review is a pretty objective proof of that. At present, it is evident that the 
key features of the underlying physics are well established and understood. At the 
same time, it is also evident that there are still many theoretical questions about 
the applications of magnetic catalysis under various conditions, where factors other 
than the magnetic field may also play a substantial role. 
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One prominent example is the dynamics of chiral symmetry breaking in QCD 
in a magnetic field. Because of a poorly understood interplay between the dynam- 
ics responsible for the quark (de-)confinement on the one hand and the magnetic 
catalysis on the other, there are a lot of uncertainties about the precise role of the 
magnetic field in this case n21lll44l l2ll25lfl46ll67ll68lR^l77^ 
[72] [34j [141 |T6l [23] |T5] |7T] . One can even suggest that there exists an intermediate 
regime in QCD, starting at magnetic fields of order B ~ 10 19 G or so, in which the 
magnetic field is sufficiently strong to provide a gluon screening 1761 and, thus, sup- 
press the vacuum chiral condensate, but still is not strong enough to produce equally 
large quark masses through magnetic catalysis 01441 . At finite temperature, further 
complications could appear because of the interplay of the magnetic field and the 
temperature in gluon screening (T6). All in all, it is obvious that there are many 
research directions remaining to be pursued in the future. 

As we argued in Sec. 13.31 magnetic catalysis may play a profound role in the 
quantum Hall effect in monolayer graphene. It appears, however, that an interesting 
variation of magnetic catalysis can be also realized in bilayer graphene [93 , 92l[86l 
l90l . In essence, it is a nonrelativistic analog of the magnetic catalysis. This fact 
alone is of interest because of a large diversity of solid state physics systems and the 
relative ease of their studies in table-top experiments. 

Finally, one should keep in mind that the fundamental studies of gauge field the- 
ories, which are known to have an extremely rich and complicated dynamics, is of 
general interest even in the regimes that are not readily accessible in current exper- 
iments. Such studies usually provide invaluable information about the complicated 
theories in the regimes that are under theoretical control. This often allows one to 
understand better the structure of the theory and even predict its testable limitations. 
In the case of QCD in a magnetic field, e.g., we may gain not only a better under- 
stand of the fundamental properties, but also get an insight into the physics in the 
Early Universe and in heavy ion collisions. 
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Appendix: Fermion propagator in a magnetic field 

Let us start from the discussion of the Dirac fermion propagator in a magnetic field 
in 3 + 1 dimensions. It is formally defined by the following expression: 

S(x,x') = i [ifd, - (JC ± ■ Y ± ) - 7C 3 f - m]' 1 <S 4 (x-x'), (80) 

where x = (x° , x 1 , x 2 , :t 3 ) = (t , r ) . By definition, the spatial components of the canon- 
ical momenta are k' = —id; — eA\ where / = 1,2,3. (The perpendicular components 
are i = 1,2.) Here we assume that e is the fermion electric charge (i.e., one should 
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take e < in the case of the electron) and use the Landau gauge A = (Q,Bx l ,Q), 
where B is the magnetic field pointing in the x 3 -direction. By definition, the com- 
ponents of the usual three-dimensional vectors A (vector potential) and r (position 
vector) are identified with the contravariant components A 1 andx', respectively. 

In the Landau gauge used, it is convenient to perform a Fourier transform in the 
time (f — t 1 ) and the longitudinal (x 3 — x' 3 ) coordinates. Then, we obtain 

S(co,p 3 ;r ± ,r' ± ) = i [y°<o - (*l ■ Y±) ~ j'p^-m] _1 <S 2 (r ± -r' ± ) 
= i[fco-(7t ± -y ± )-y 3 p 3 + m\ 

x [(0 2 -1t\ + ieBYY-(p 3 ) 2 -m 2 }~ 1 S 2 (r ± -r' ± ), (81) 

where rx is the position vector in the plane perpendicular to the magnetic field. 

In order to obtain a Landau level representation for the propagator dSTJ, it is 
convenient to utilize the complete set of eigenstates of the operator lt\. This oper- 
ator has the eigenvalues (2k+ l)\eB\, where k = 0,1,2,... is the quantum number 
associated with the orbital motion in the perpendicular plane. The corresponding 
normalized wave functions read 



Wk P2 (r±) 



1 



1 



-Mk 



X j+P2e)e-i 



-is ± .x P2 



(82) 



where H^z) are the Hermite polynomials Ii94| , I = 1 j^J\eB\ is the magnetic length, 
and s± = sign(eB). The wave functions satisfy the conditions of normalizability and 
completeness, 



d 2 r±W*k P2 {r±)¥k' P ' 2 ^±) = &kk>8{P2- p' 2 ), 



(83) 



dpi £ Wk P2 (r±)Wk P2 ( r '±) = 5 2 (r±-r' ± ), 



k=0 



(84) 



respectively. 

By making use of the spectral expansion of the 5-function in Eq. (1841) . as well as 
the following identities: 



(*±-r±)Vkpi = ~ ( f I V 2 (* + 1 ) Wk+ 1 , P1 9- - V2k }jr k - 1 <P2 



2k+l 

p. Wkp 2 



(85) 
(86) 



with !?± = j (l ± is±y l y 2 ) being the spin projectors onto the direction of the mag- 
netic field, we can rewrite the propagator in Eq. ( TSTb as follows: 



S(a>,p 3 ;r ± y ± ) = dp 2 £ i[fco- (* x -yJ - yY + m] co 2 -( P 3 ) 2 



k=0 
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-(2k+ l)\eB\ +ieB Y Y -m 2 ] ~ V 2 ( r ±M* P >x) 
= e i ^ r ^ r '± ) S(co, P 3 ;r ± ~r' ± ). (87) 
The Schwinger phase is given by 

( X l +x n ){x 2 -x l2 ) 

and the translationary invariant part of the propagator reads 

y ± .(r ± -r' ± )L, 1 ,_ 1 (^), (90) 



where we used the short-hand notation 

q 2l 2 ' K ' 

In order to integrate over the quantum number p2 in Eq. (|87| |. we took into account 
the following table integral J94): 

e- x2 // m (x +y)/f„(x + = 2"K^ 2 mlz"- m L'^ m (-2yz), (92) 



which is valid when m<n. Here L" are the generalized Laguerre polynomials, and 

— ■ 

Here a short remark is in order regarding the general structure of the Dirac prop- 
agator in a magnetic field. It is not a translationally invariant function, but has the 
form of a product of the Schwinger phase factor e ! *( r -L> r jJ and a translationally in- 
variant part. The Schwinger phase spoils the translational invariance. From a physics 
viewpoint, this reflects a simple fact that the fermion momenta in the two spatial di- 
rections perpendicular to the field are not conserved quantum numbers. 

The Fourier transform of the translationary invariant part of the propagator d89l 
reads 

Zfr, 3 . _ \ „--p^f- (-l)"A,(fi),p 3 ;px) 

S(co,p ,p ± )=2* e p± Z—^-——, (93) 



where 



D„(o),p 3 ;p ± ) = (fco - rV +m) [L„ {2p\i 2 ) & + -L„_i (2/^ 2 ) 

+2(Y ± -p ± )L l n _ l (2pie 2 ). (94) 



Magnetic Catalysis: A Review 



29 



Taking into account the earlier comment that the perpendicular momenta of charged 
particles are not conserved quantum numbers, this representation may appear sur- 
prising. However, one should keep in mind that the result in Eq. (|93l l is not a usual 
momentum representation of the propagator, but the Fourier transform of its trans- 
lationary invariant part only. 

In some applications, it is convenient to make use of the so-called proper-time 
representation 111661 . in which the sum over Landau levels is traded for a proper-time 
integration. This is easily derived from d93l i by making the following substitution: 



(95) 



C0 2 -2n\eB\ - (p 3 ) z -m 2 + i0 



= I d ie "K-2»|eB|-(p 3 ) 2 -m 2 +i-0] 



Then, the sum over Landau levels can be easily performed with the help of the 
summation formula for Laguerre polynomials |94| . 



£L°(*)z" = (l- Z )-( a+1 )exp 



n=0 



XZ 

z-l 



(96) 



The final expression for the propagator in the proper-time representation reads 



5(c0,/? 3 ;p ± ) 



dse 



is[a> 2 -m 2 -(p 3 ) 2 ]~i(p 2 ± e 2 )tan(s\eB\) 



y°co-(y-p)+) 



"(pV - p 2 Y l )tan{seB) [l - y 1 y 1 tan(seB)\ , 



(97) 



where (y-p) = (y ± -p±) + yV- 

Using the same method, one can also derive the Dirac fermion propagator in a 
magnetic field in 2 + 1 dimensions. It has the same structure as the propagator in 
Eqs. d87]i, <[88]>, (89), and <|90]>, but with p 3 = 0, i.e., 



where 



5 2+ iO;r-r') = i 



5 2+I (co;r,r')= e ,<i,(ry) 52 + i(co;r-r'), 
fco + m 



(98) 



2k£ 2 \ 
i y-(r-r') 



CO 2 — 2n\eB\ -m 2 



L„(^)^-+L„_i(^)^+] 



(99) 



i 2 (0 2 -2n\eB\-m 
The Fourier transform of the translationally invariant part is 

\f(0 + m) [L n (2p 2 £ 2 ) -L„_] (2p 2 i 2 ) &J\ 



(y-p) 



(D 2 -2n\eB\ 



(O 2 -2n\eB\ —m 



2*2\ 



(100) 
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Finally, the proper-time representation reads 



5 2+ i(<o;p) = / dse 



!> is[(0 2 -m 2 }-i(p 2 e 1 )tnn(s\eB\) 



y° co- (y-p)+j 



(p x y 1 -p 1 Y l )iw(seB) [l - yV tan(seB)] . (101) 
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